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diffusion, convection, and adiabatic deceleration under the assumption of a spherically symmetric 
heliosphere. The model utilizes a comprehensive database of GCR measurements from various 
particle detectors to determine boundary conditions. By using an updated GCR database and 
improved model fit parameters, the new BON model (BON14) is significantly improved over the 
previous BON models for describing the GCR radiation environment of interest to human space 
flight. 
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1. Introduction 


It is well known with many years of research and as documented in the literature that the ioniz- 
ing nature of Galactic Cosmic Ray (GCR) particles poses a potential health risk for crew members 
in space, particularly for future long-term missions in free-space [1-3]. Long-term exposure to 
mixed (both low and high) Linear Energy Transfer (LET) GCR radiation significantly increases 
the Risk of Exposure Induced Cancer (REIC) and Risk of Exposure Induced Death (REID) [1-3]. 
Another significant concern arises from the interaction of GCR particles with the electronics inside 
and outside of a spacecraft. Energetic GCR particles may deposit energy in electronics, e.g., micro- 
processors, memory units, sufficient to cause memory bit flips and latch-up, which arc generically 
called Single Event Effects (SEE) [4,5], 

In order to evaluate the potential risks induced by the GCR ions, the Badhwar-O' Neill (BON) 
GCR flux simulation model [6-10] has been developed to numerically solve the Fokker-Planck 
(FP) equation. The BON model takes into account diffusion, convection, and adiabatic deceleration 
within the heliosphere and provides the flux of GCR particles of a given charge, Z, as a function 
of energy near earth ~ 1 astronomical units (AU) in free-space beyond the Earth’s magnetosphere. 
The solution is obtained under the assumptions of a quasi-steady state and a spherically symmetric 
interplanetary medium [11, 12]. With these assumptions, the FP equation can be written as [13]: 
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where r is the radial position in units of AU; T is the kinetic energy (MeV/n); U(r, T) is the GCR 
flux, V s (r) the solar wind speed (~ 400 km/s); k(i\ T ) the particle diffusion coefficient tensor; and 
ct(T) = (T + 2Eo)/(T + £b), with Eq being the rest energy per nucleon (E p ~ 938 MeV/n) of the 
GCR particle. The solution also assumes that at a boundary distance r = R/,. modulation of U (r, T) 
is negligible, and therefore provides the boundary condition at U ( R /, . T ) = l Jo as a known quantity. 
This quantity, Uq, is ion specific and parametrically described by several free parameters, which 
arc known as Local Interstellar (LIS) parameters. 


1.1 Description of the LIS parameters in BON14 

At a distance well outside of the solar system, around R/, = 100 AU, the GCR modulation due 
to the turbulent solar wind and heliospheric magnetic field is negligible and therefore each GCR 
ion energy-flux spectrum is assumed to be constant. This constant GCR field is referred to as the 
LIS flux spectrum (Uq) and represents one of the boundary conditions for the BON. The LIS flux 
for each GCR ion has the following parametric relationship as a function of charge (Z) and kinetic 
energy: 


Uo(Z,T) \ Rb=mAU =jo(Z)(T N + E 0 y {z) p^ 1 p S ^ z \T + E 0 )~ r{z] (1.2) 

where jo, 8 , and y arc free parameters for each GCR ion, /3 = v/c is the velocity of the ion relative 
to the speed of the light, and /jy is the relative velocity at 7y = 35 GeV/n (selected arbitrarily for fit- 
ting). 
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The LIS parameters are formulated by using the GCR mea- 
surement data from detectors at or near 1 AU, e.g., satellite 
and balloon measurements. In the model, the flux of any 
ion beyond nickel (Z > 28) is obtained by scaling from the 
silicon result. In BON14, as opposed to the previous BON 
releases, we have modified the LIS parameters, jo, 8, and y, 
driven by a sensitivity study by using several metrics. For 
a detailed description of the sensitivity analysis, please see 
the referenced publication series by T. Slaba et al. [14-16], 

One of the major results of this study was that for the 
differential effective dose rate as a function of kinetic en- 
ergy behind 20 g/cm 2 of aluminum shielding at a period of 
minimal sun activity, GCR ions in the energy domain be- 
tween 0.5 GeV/n and 4.0 GeV/n account for most of the 
exposure. As presented in the study, H and He with bound- 
ary energy (BE) less than 0.5 GeV/n induce approximately 
9% of the total effective dose. GCR ions with Z > 2 and 
BE less than 0.5 GeV/n induce less than 4%, while all GCR 
ions, from hydrogen (Z=l) to nickel (Z=28), with BE be- 
tween 0.5 GeV/n and 4.0 GeV/n induce ~ 66% of the total 
effective dose behind 20 g/cm 2 aluminium shielding. 

Based on the results of that study, the new LIS parame- 
ters are fitted to the GCR data, in such a way that allowed us 
a range of parameter combinations to be evaluated to min- 
imize the uncertainty and relative difference between GCR 
measurements and the model. The updated LIS parameters 
for each of the corresponding ions from hydrogen to nickel 
are shown in Table 1 . 

1.2 Selection of GCR data 

In this revision of the BON model, BON14, we have 
included the GCR data beyond 1970, which spans Solar 
Cycles 20 to 24 (to date). A comprehensive table that rep- 
resents the entire GCR measurements used for the BON 14 


Table 1: The LIS parameters that are 
used in the BON14 shown for each 
ion. 


z 

7 

8 

jo 

1 

2.75 

-2.82 

9.50xl0~ 4 

2 

2.80 

-2.00 

4.53 xl0~ 5 

3 

3.21 

-0.69 

6.37 xl0~ 8 

4 

2.93 

1.50 

1.20xl0~ 7 

5 

3.00 

-0.40 

2.40x 10 -7 

6 

2.70 

-2.00 

1.60xl0~ 6 

7 

2.95 

-0.60 

2.65 xlO -7 

8 

2.73 

-1.90 

1.50xl0~ 6 

9 

3.08 

0.40 

1.63xl0~ 8 

10 

2.75 

-1.60 

2.35 x 10 -7 

11 

2.73 

-1.80 

4.60xl0~ 8 

12 

2.70 

-2.40 

3.03 xlO -7 

13 

2.75 

-1.40 

5.30xl0~ 8 

14 

2.65 

-2.40 

2.65 x 10 -7 

15 

3.15 

2.00 

5.68xl0~ 9 

16 

2.70 

-1.00 

5.78xl0~ 8 

17 

3.13 

2.00 

5.99xl0~ 9 

18 

2.90 

0.60 

1.68xl0~ 8 

19 

3.13 

0.80 

7.90xl0~ 9 

20 

2.75 

-1.60 

3.23 xlO -8 

21 

3.15 

0.40 

3.50xl0~ 9 

22 

3.00 

-0.50 

1.44 xl0~ 8 

23 

3.00 

-0.50 

7.14x 10~ 9 

24 

2.90 

-1.00 

1.78xl0~ 8 

25 

2.80 

-1.00 

1.39xl0 -8 

26 

2.60 

-2.40 

2.00xl0 -7 

27 

2.60 

-2.50 

l.llxlO -9 

28 

2.55 

-2.40 

1.19xl0~ 8 


is presented in [10, 16]. In the past, the LIS parameters of 

the BON were uniquely influenced by measurements from the Cosmic Ray Isotope Spectrometer 
(CRIS) on the NASA Advanced Composition Explorer (ACE) spacecraft [17]. The CRIS instru- 
ment is currently measuring the flux of ions and their isotopes from boron (Z=5) to nickel (Z=28), 
where the lowest and highest kinetic energy measurements are ion specific. CRIS provides ki- 
netic energy of GCR isotopes between ~ 50 - 500 MeV/n. In the BON14, greater emphasis was 
placed on the higher kinetic energies, a region not covered by CRIS. Nonetheless, we show that the 
updated model still accurately represents the lower energy regions. 
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2. Comparison of BON14 to GCR measurements 

In Fig. 1, the differential flux for H ion as a function of energy for several GCR measurements, 
AMS (Alcaraz et al„ 2000) [18], IMAX (Menn et al„ 2000) [19] , and PAMELA (Adriani et 
al., 2011) [20], compared with the BON14 model for the same measurement periods is shown. 
The solid black line in the figure (LIS-BON14) represents the LIS flux for H. The energy axis is 
magnified in the 0. 1- 10 GeV/n region and re -plotted in the inset figure. 



Figure 1: The differential flux for H ions as a function of energy is shown for various GCR measurements, 
AMS (1998), IMAX (1992), and PAMELA (2008). Those measurements are compared with BON14 for the 
same measurement periods. The solid black line (LIS-BON14) represents the LIS flux for H ion at R=100 
AU. The energy axis is magnified in the 0.1-10 GeV/n region and replotted in the inset figure. 


In this paper, we report the comparison of BON 14 with the GCR measurement data by evaluat- 
ing the relative differences between both measurements and the model. A comprehensive study in- 
cluding model uncertainty analysis comparison of the previous BON models (BON10 and BON1 1) 
with the current model was reported elsewhere [10]. The relative difference, Rd, is a measurement 
between the GCR measurement data and the BON model that allows us to find where the model 
underpredicts or overpredicts the measured GCR flux. Rcl is defined as follows: 


1 N 

Rd = V 

k= l 


Models — Deltas 
Datcik 


(2.1) 


where N is the number of measurements. The average of the residual Rd does not provide a 
complete picture of the model versus measurement agreement, as the quantity may end up in the 
vicinity of zero due to positive-negative cancellation effect. Therefore, we also used the average 
absolute relative difference, \Rd\, to determine the overall difference between the model and the 
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Figure 2: The 27-day average differential flux for O ions at 230.8 MeV/n from ACE/CRIS (circle) as a 
function of date is compared with BON14 (solid line). On the right axis their Rd is shown (diamond). 


GCR measurement database, which can be described as follows: 


|JM| 



| Mod elk — Datcik \ 
Data b 


( 2 . 2 ) 


The latter metric better quantifies the spread in model errors and may be used as a quality check for 
the model. In Fig. 2, the differential flux averaged over 27 days for O ions at 230.8 MeV/n from 
ACE/CRIS (circle) as a function of date is compared with BON14. For this particular energy bin, 
\Rd\ is calculated as 13.6%. 

In addition, we compared BON14 with the previously released version, BON1 1, in the energy 
region of interest. In Fig. 3, comparisons of the average of the integrated differential flux for H, 
Fie, O, and Fe ions for all GCR data, excluding ACE/CRIS data, with BON11 and BON14 models 
are shown. In these figures, we separated the data into two energy regions to simplify comparisons, 
where they represent lower energy bins (< 4.0 GeV/n) and higher energy bins (> 4.0 GeV/n). In 
those plots, we define the average flux (< U >) as the integral of the differential flux (dU /dE) over 
all energy bins, which is given by the following equation: 


1 fdU IT , 
< C >= — — / — dE, 


N b , 


dE 


(2.3) 


where N b - in is the number of energy bins and dE is the bin width of the data set. Average flux was 
evaluated by summing over the reported energy bins. The \Rd\ values for low-energy H ion are 
18.7% and 13.6% for BON11 and BON14, respectively. Moreover, as shown in the Fig. 3 for Fie, 
O, and Fe, BON14 presents a significant improvement over BON1 1 especially for energies below 
4.0 GeV/n. 
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In Table 2, the overall average relative difference metrics, Rd and Rd | , are presented for 
BON11 and BON14. These values were determined by averaging the relative difference metrics 
over each ion. Individual results for each ion were reported in [10]. The results in Table 2 were 
obtained by considering the entire GCR measurement database (all ion and energy bins). It is seen 
in Table 2, that for BON11, the \Rd\ value is calculated as 23.7% and Rd = 17.9%. For BON14, 
these errors are reduced to \Rd\ = 13.0% and Rd = —0.4%. The BON11 model systematically 
overpredicts the GCR measurement data, whereas BON 14 provides a more balanced prediction. 
More importantly, the overall spread in the new model error was reduced from 23.7% to 13.0%. It 
should be emphasized that at all energies, BON 14 is only a marginally improved fit to the hydrogen 
GCR data compared with BON11. There is room for improvement since BON14 underpredicts 
hydrogen GCR data for ~ 8 year-long periods since 1998 (per Fig. 3). A similar underpredicting 
pattern is observed for He as well. 

3. Conclusions 

We present the comparison results of BON14 with an 
updated set of GCR measurements from various balloon, 
satellite, and space shuttle measurements. The new model 
LIS parameters were fitted with an updated approach based 
on the sensitivity study described earlier. This study showed 
that the GCR ions with energies between 0.5 GeV/n and 4.0 
GeV/n account for most of the shielded effective dose. 

The overall average Rd value is calculated as 17.9% 
for the BON11 and -0.4% for the BON14 model, whereas 
\Rd\ is found as 23.7% and 13.0% for BON1 1 and BON14 
models, respectively. As described in [10], the uncertainty metrics study also revealed that overall 
the BON14 model has been improved significantly with respect to the previous models (BON10 
and BON11) both at low- and high-energy regions. However, it should be noted that the LIS 
parameters used by BON14 slightly undeipredict hydrogen GCR data for most of the years since ~ 
1998 (beginning of a solar m a ximum). We anticipate significant model improvement in the energy 
region of interest with hydrogen data when the new GCR measurements are available in the future, 
e.g., AMS-02 measurements. In addition, we will attempt to improve the model by incorporating 
several physics effects on the GCR particles (e.g., curvature and gradient drifts) which have not 
been included in the previous BON models. 


Table 2: The average relative differ- 
ence metrics between the entire GCR 
database with BON11 and BON14. 



Average [%] 


Rd \Rd\ 

BON11 

17.9 23.7 

BON14 

-0.4 13.0 
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Year Year 


Figure 3: The average flux values, < U >, for H ion: (a) below, (b) above; He: (c) below, (d) above; O: (e) 
below, (f) above; and Fe: (g) below, (h) above 4.0 GeV/n are shown as a function of measurement date for 
BON11 and BON14. |/?c/| values of the models compared against same GCR database are presented in the 
legends of the figures. 
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